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Abstract 


The paper gives the theory of magnetic propulsion of liquid lithium streams 
and their stability in tokamaks. In the approximation of a thin flowing layer 
the MHD equations are reduced to one integro-differential equation which 
takes into account the propulsion effect, viscosity and the drag force due to 
magnetic pumping and other interactions with the magnetic field. A criterion 
is obtained for stabilization of the “sausage” instability of the streams by 


centrifugal force. 
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For fusion, liquid lithium has unique physical properties (low melting temperature of 
Tmett ~ 180° C, low density p œ 0.53 g/cm?, low viscosity v ~ 0.0004 Pa- sec in combination 
with high thermal capacity cp ~ 4200J/ (kg: K) and thermal conductivity x ~ 47 W/(m- K) 
[1]) which would make it potentially the best coolant for blankets of thermonuclear devices. 
Nevertheless, because of the rather high electric conductivity, o ~ 3.3. 106 (Q . m) t [2]) 
liquid lithium has been essentially dropped from the design of magnetic fusion reactors. In 
the strong magnetic field of tokamaks it is very energy consuming to drive a lithium flow 
if it is restrained by pipe walls or channels. Existing proposals for using so-called “liquid 
walls” for power extraction from the tokamak reactor plasma [3] rely on passive lithium flow 
which is sensitive to a variety of MHD effects. A recently discovered (December, 1998) effect 
of magnetic propulsion suggests another way of using lithium, i.e., with the actively driven 
free surface streaming along the inner walls of the vacuum chamber. 

This paper outlines the theoretical basis of magnetic propulsion. Fig. la shows two 
lithium layers of about 1 cm thick, which are held to the surface by the Ipo X Btor electro- 
magnetic force, where Ipa represents an externally driven electric current through the layers 
and Bio, is the toroidal magnetic field of the tokamak. The pressure p in the fluid, created 
by the electromagnetic force, is inhomogeneous in both normal and longitudinal directions. 
In the normal direction it is analogous to the atmospheric (highly magnified) pressure of 
the Earth. It vanishes at the free surface of the liquid and reaches maximum (about one 
or several atmospheres) at the surface of the wall. Its inhomogeneity along the wall surface 
propels the fluid (as a wind) from the high field side of the chamber, where the Ipo; x Bior 
force is largest, to the low field side. 

Neglecting viscosity and possible electromagnetic drag forces, the fluid velocity V in each 
flowing layer obeys the Bernoulli law 

y2 
p + pgz + Pog = const, p= |Lw xX Biorl, (1) 
where z is the vertical coordinate and g = 9.8 m/sec? is the acceleration of gravity. Note, 


that the tokamak magnetic field Bir = 5 T, reasonable for reactors, corresponds to a 
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magnetic pressure of 100 atm. Using only a small fraction of it (e.g., 1-2 %) it is possible to 
create a pressure drop along the flow larger than 1 atm, which corresponds to ~ 20 m/sec 
of the liquid lithium velocity. Such a velocity is sufficient to provide reactor relevant power 
extraction from the plasma of the order of 3.5 MW /m?. For comparison, in the International 
Thermonuclear Experimental Reactor project this number is about 20 times smaller [4]. 

This paper describes the propulsion itself. Two closely related problems of how to supply 
the liquid to the inlet (e.g., in the form of fast thin jets and their conversion into the flow) 
or to remove it from the chamber (e.g., by converting the flow into a jet shower) are not 
considered here. 


The dynamics of the liquid metal is determined by a standard set of MHD equations 


V A 
Pr +p(V-V)V = -Vp — pgVz + (j x B) + vAV, (2) 
ðA j 
ial ee B)= 2 


where A is the vector potential of the magnetic field, yx is the electric potential, j is the 
current density in the metal and v is the viscosity. 

The toroidal geometry of tokamaks justifies the consideration of axisymmetric flow with 
no side walls limiting the streams. Note, that even with a few separating walls present, the 
magnetic propulsion is not sensitive to the resulting Hartmann drag force [5]. 


We use the orthogonal curvilinear coordinates n, s, ó 
/ 
T = Tuy — Zn, Z= Zw ae (4) 


suitable for describing the thin liquid flow along a guide plate as is shown on Fig. 1b. 
Here, r, ó, z are cylindrical coordinates, the functions rw = rw(s), zw = Zw(s) describe the 
geometry of the wall surface in terms of the length s along the surface, r’? + 42 = 1. The 
normal coordinate n is 0 at the guide surface and n = h at the free surface of the flow 
(here, h = h(s) is the thickness of the metal layer). It is also convenient to introduce the 
normalized coordinate 7 (conformal to the flow surface) and the poloidal curvature & of 


the guide wall surface 


ri 
ww — w wl 


N=—-, &Kew& |r 
In thin layer approximation we adopt the relations 
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We also can neglect the normal component of velocity 
V=VVs, VVn & V. 
The height h of the stream is determined by the continuity equation 
‘ h 1 
M = anp f rVdn = 2rph | TV án, 
0 0 


where M is the mass flow in the stream. 
The electric current can be represented through its stream functions as 


; ! F; ! F; 
HoJ = (VEF x V¢), Jn = > Js = 
HoT 


Ugh’ 


The stream function of the electric current F includes three components 


F = Fy + Fenn + iln, s), Fy =T Bior = const. 


(10) 


The major term F,, = 0.2/7rc is related to the total current Jprc in the toroidal field 


coils, Feztñ represents the poloidal current (uniformly distributed across the stream), Fert = 


0.27, driven through the stream by an external source and 7 is the stream function of the 


eddy currents. Assuming wall thickness smaller than h and neglecting the current sharing 


between the wall and the stream, i satisfies the boundary conditions i(0, s) = i(1, s) = 0. 


In thin layer approximation the MHD equations can be written in terms of their projec- 


tions as 
OV, ph $ FF! 
= wV? = — — z ; 
Fear, He h Oh hin? 
OV FF' v 
Pap + PVV; = =p, — PIžw — E + a Van; 
l, VF i F; 


mana. 


(11) 
(12) 


All s-derivatives here are taken at n = const. Eqs. (11-13) together with Eq. (8) represent a 
full set of MHD equations. We kept the time derivative of the normal component of velocity 
in the first equation in order to cover also possible instabilities of the streams. Assuming 
that the flow is directed along the magnetic surfaces of the equilibrium configuration, we 
neglected here the interaction of the flow with the poloidal magnetic field. It can be easily 
included as an additive effects. 

The normal component of Ohm’s law in the leading approximation gives the distribution 


of the scalar potential 


UF y n k 
PE = pols) — ae} =- hrVdñ, (14) 
while the longitudinal component 
UN” Fenti 
I F, G) = ext ñ 15 
Po Pos Loo hr ey 


determines the stream function (7, s) and eddy currents 
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. 1 N ~ sy 1 
. =-Fu(=/ Can) =v- | Uda, (16) 
Loohr r 0 s 0 
F, 1 m al U i 
r r Jo sls POPS 


The normal component of the equation of motion in leading order (neglecting the cen- 


trifugal force) gives the pressure distribution inside the fluid 


F?(1) — F? Feal- n) —$ 
- (1) =F, (1 — 7) fi. (18) 
2H0T” Hor? 


Its substitution into the longitudinal equation of motion leads to a magnetic propulsion 


equation (MPE) 


> 
+ pVVi == p97, = 


POE [rpcuo r? 
1 Pa sae on VA v 
2 = 
—ghr Fá (= if Wan) (=). + Vin 
—(15B;),(1 = nr) = oV B?, (19) 


which contains only one unknown variable V (7, s). 
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'The first term in the right hand side of MPE is gravity, the second is the magnetic 
propulsion term, which accelerates the fluid toward the low magnetic field side of the guide 
wall. The third term in the right hand side of MPE represents the diamagnetic drag force or 
“magnetic pumping”, the fourth is viscosity. The last two terms, written without derivation 
describe interaction with the poloidal magnetic field. The first of them is the propulsion 
effect due to interaction of a toroidal current Jy in the stream, if present, with tangential 
component of the poloidal field B,, while the second is the drag force due to interaction 
with the normal (axisymmetric) B, component of the poloidal magnetic field (e.g., due to 
plasma displacement). 

The diamagnetic drag force with respect to the propulsion term can be assessed assuming 
a uniform velocity profile across the stream. For values averaged over cross section, their 


ratio is determined by the propulsion parameter Rp as 


Irre 
~ Rp=R 20 
propulsion a. 2 i ma (20) 


pga h?V 


drag force Ror, 


Rə (21) 


The same propulsion parameter is responsible for possible peeling of the flow surface due to 
the appearance of “negative” pressure in the flow, if at some point Rpr’, > 1. 

For liquid lithium poo ~ 4 and for typical rector parameters, e.g., r = 6 m, V œ 20 
m/sec, h ~ 1 cm, the diamagnetic drag force has only a small effect on flow. Even for 
Ipa = 0.01Irrc, the propulsion parameter Rp < 1. It is important to notice that the 
Reynolds number Rə, introduced here, is different from the “electromagnetic” Reynolds 
number Rı = uocghV, which plays a role in plasma stabilization by the lithium streams 
[6,7]. This fact makes lithium dynamics consistent with the possible plasma stabilization by 
intense lithium streams. 

The interaction with the normal component of the magnetic field is characterized by the 


conventional Reynolds number Ro 


—- ds, Ro z= Loo LV, (22) 
Ho 


where L is the length of the flow. Because Ro >> 1, the flow is very sensitive to the quality of 
alignment of the flow with magnetic surfaces. Nevertheless, the magnetic propulsion allows 
a certain level of misalignment of the flow with the magnetic surfaces and of the normal 
magnetic field. 

The possible electromagnetic instabilities appear in the next approximation when the 
centrifugal force and the dependence of the magnetic field line tension on the major radius 


are taken into account in the equation for pressure 


, 
(» — oS"). = —pgrh — Pah — PkyhV? (23) 
due to explicit dependence of the RHS on the height h of the flow. Thus, the first term 
drives the gravitational instability (at the top of the vacuum chamber) [8], while the second 
one drives the “sausage” instability at the high field side, where &, < 0 [9]. The third term 


is the stabilizing effect of the centrifugal force. 


The centrifugal force robustly stabilizes the sausage instability of the stream if 


nr 1 y2 i 
| PKu V? — 2p-. | dn > 0, i > ay. (24) 
-1 r 2 Q2KyT 


where py is the pressure of the stream on the wall surface. Because |z/,| < 1 and the 
poloidal curvature is determined by the minor radius a of the plasma, kw œ 1/a, this 
stability condition can be fulfilled for sufficiently high inlet stream velocity. 

Much smaller flow speed is required for stabilization of the gravitational interchange 


instability 


(v?) > a (25) 


for the stream moving along the upper half of the vacuum chamber. 

An example of a stationary solution to MPE (19) for liquid lithium is shown on Fig. 2. 
The velocity, averaged over the cross-section, and shown on Fig. 2a, rises from an inlet value 
to the characteristic level of 20 m/sec at the middle of the path. The dotted line solution 
was obtained neglecting viscosity and the diamagnetic drag force in MPE. Neglecting dia- 


magnetic drag force (dot-dashed line) shows the very small effect of diamagnetic drag (in 
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contrast to the case of Ref. [10], where this effect was exaggerated). The normal component 
of magnetic field does affect the propulsion, as is shown by a dashed line solution. Note, 
that in this case the imposed vertical magnetic field would correspond 2-3 h(0) in the plasma 
displacement (depending on the plasma current). 

The flow becomes significantly thinner on the way to the outlet point (Fig. 2b), Fig. 2c 
shows a comparison of the velocity profile for the cases without the effect of diamagnetic drag 
force, without both viscosity and drag force and with a normal magnetic field. Fig. 2d shows 
a comparison of the driven current density (solid lines) with the eddy currents (dot-dashed 
lines), which are significantly smaller at all cross-sections of the flow. 

The effect of viscosity on magnetic propulsion in the laminar approximation is very weak 
and for typical velocities 10-20 m/sec of the streams is localized in a sub-millimeter boundary 
layer near the surface of the wall. Note that, while magnetic propulsion drives the inner 
layer of the flow, viscosity propels the surface layer. 

Summarizing, the magnetic propulsion of intense lithium streams represents an means 
for driving controlled lithium streams in a tokamak magnetic field. With development of the 
technology of injection/ejection of the streams it can be used for efficient power extraction 
from the plasma. The magnetic propulsion equation (19) describes the 2-dimensional MHD 
of the liquid lithium flow and shows robustness of the propulsion effect to viscous drag and 
diamagnetic drag force, characterized by dimensionless parameters Rp < 1 and R < 1. 
While the flow is sensitive to the normal components of the magnetic fields, in the case of 
magnetic propulsion misalignment by a few thicknesses of the flow is still tolerable. We also 
have shown that the streams can be stabilized by centrifugal force against “sausage” and 
gravitational instabilities. The thin layer approach can be extended to cases which would 
include disturbances of axisymmetry (like penetrations or non-symmetric normal fields). 

This work was supported by United States Department of Energy Contracts No. DE- 
AC02-76-CHO-3073. 
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FIGURES 
FIG. 1. (a) Cross-section of a tokamak with the Toroidal Field Coils (TFC), a vacuum chamber 


with two lithium streams at the top and bottom halves of the vacuum chamber. (b) local coordinate 


system n, s. 


FIG. 2. Numerical soluton of the MPE for R = 6 m, a = 1.6 m, Bitor = 5 T (at r = R), 
Iet = 0.01Irrce = 1.5 MA, hiniee = 2 cm. (a) averaged velocity and (b) height h as function of 
longitudinal coordinate, (c) flow velocity profile in 3 cross sections along the stream. Solid lines 
represent solutions of MPE (19) for B, = 0, dashed lines for Bn = 0.005Bjor cos(s/a), dotted line 
corresponds the “ideal” solution with neither viscosity nor magnetic pumping, and dot-dashed lines 
the solution with no magnetic pumping. (d) profiles of externally driven current density (solid line) 


and of eddy current (dashed line). 


10 


axis of symmetry 


J x B Force 


V, 


stream 


Li stream 1 


11 


<V>, m/sec 


12 


0 
0 2 6 8 s/L 1 
eu . 1h i 
Inlet 8 1/2 path g ba Sg, Outlet 
6 6 Sg 
4 4 
.24 .24 
| , : ' 0 - zes Lit 0 : i i 
0 10 20 30 0 10 20 30 0 10 20 30 
V, m/sec V, m/sec V, m/sec 
n/h n/h n/h 
10 a 12% 
a| 14 path 8 i 1/2 path 8 i 3/4 path 
6 4 6} | 6 | 
Ah Ah Hi 
2 | 2} fi 2} f 
o+ a ol ol | 
0 4 6 0 2 4 0 2 4 
j, MA/m*2 j, MA/m^2 j, MA/m^2 


External Distribution 


Plasma Research L aboratory, Australian National University, Australia 
Professor I.R. J ones, Flinders University, Australia 

Professor J oao Canalle, Instituto de Fisica DEQ/IF - UERJ, Brazil 
Mr. Gerson O. Ludwig, Instituto Nacional de Pesquisas, Brazil 
Dr. P.H. Sakanaka, Instituto Fisica, Brazil 

The Librarian, Culham Laboratory, England 

Library, R61, Rutherford Appleton Laboratory, England 

Mrs. S.A. Hutchinson, J ET Library, England 

Professor M.N. Bussac, Ecole Polytechnique, France 

Librarian, Max-Planck-Institut fur Plasmaphysik, Germany 

J olan Moldvai, Reports Library, MTA KFKI-ATKI, Hungary 

Dr. P. Kaw, Institute for Plasma Research, India 

Ms. P.J. Pathak, Librarian, Insitute for Plasma Research, India 
Ms. Clelia De Palo, Associazione EURATOM-ENEA, Italy 

Dr. G. Grosso, Instituto di Fisica del Plasma, Italy 

Librarian, Naka Fusion Research Establishment, J AERI, J apan 
Library, Plasma Physics Laboratory, Kyoto University, J apan 


Research Information Center, National Institute for Fusion Science, J apan 


Dr. O. Mitarai, Kyushu Tokai University, J apan 


Library, Academia Sinica, Institute of Plasma Physics, People’s Republic of China 


Shih-Tung Tsai, Institute of Physics, Chinese Academy of Sciences, People’s Republic of China 


Dr. S. Mirnov, TRINITI, Troitsk, Russian Federation, Russia 
Dr. V.S. Strelkov, Kurchatov Institute, Russian Federation, Russia 


Professor Peter Lukac, Katedra Fyziky Plazmy MFF UK, Mlynska dolina F-2, Komenskeho 


Univerzita, SK-842 15 Bratislava, Slovakia 
Dr. G.S. Lee, Korea Basic Science Institute, South Korea 
Mr. Dennis Bruggink, Fusion Library, University of Wisconsin, USA 
Institute for Plasma Research, University of Maryland, USA 
Librarian, Fusion Energy Division, Oak Ridge National Laboratory, USA 
Librarian, Institute of Fusion Studies, University of Texas, USA 


Librarian, Magnetic Fusion Program, Lawrence Livermore National Laboratory, USA 


Library, General Atomics, USA 


Plasma Physics Group, Fusion Energy Research Program, University of California at San 


Diego, USA 
Plasma Physics Library, Columbia University, USA 


Alkesh Punjabi, Center for Fusion Research and Training, Hampton University, USA 
Dr. W.M. Stacey, Fusion Research Center, Georgia Institute of Technology, USA 
Dr. J ohn Willis, U.S. Department of Energy, Office of Fusion Energy Sciences, USA 


Mr. Paul H. Wright, Indianapolis, Indiana, USA 


03/26/01 


The Princeton Plasma Physics Laboratory is operated 
by Princeton University under contract 
with the U.S. Department of Energy. 


Information Services 
Princeton Plasma Physics Laboratory 
P.O. Box 451 
Princeton, NJ 08543 


Phone: 609-243-2750 
Fax: 609-243-2751 
e-mail: pppl_info@pppl.gov 
Internet Address: http://www.pppl.gov 


